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Light Quark Physics with Dynamical Wilson Fermions* 
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We present results for spectroscopy, quark masses and decay constants obtained from SESAM's and T^L's large 
statistics simulations of QCD with two dynamical Wilson fermions. 



1. Introduction 

The generation of gauge configurations with 
two dynamical Wilson fermions, started by the 
SESAM -collaboration in spring '95 and expanded 
by the formation of the TxL -collaboration in '96, 
is now well under way. The primary goal, the 
generation of a set of configurations comprising 
C(5000) trajectories at three sea-quark masses 
with m^/mp = 0.84,0.76,0.69 at /? = 5.6 on 
16^ X 32 lattices was achieved well on time at the 
end of '96. Given the rather small box-size of this 
simulation the aim of the T^L-collaboration is to 
study finite size effects by simulating at the light- 
est of SESAM's mass- values but on a 24^ x 40 lat- 
tice. In addition, we are pushing to even smaller 
quark masses, hopefully close to m^/mp ~ 0.5. 
Last but not least, in an attempt to clarify is- 
sues surrounding the approach to the chiral limit, 
we have added a fourth sea-quark value to the 
16^ X 32 simulation and the data acquisition for 
this mass has just been completed. 

In this short note, I summarize some of the 
main findings of our analysis of two-point cor- 
relators on these gauge configurations with the 
emphasis on the methods we have applied; re- 
sults are preliminary and are to be updated soon, 
when data at the additional sea-quark mass at 
rriTr/mp = 0.81 will be included in the analysis. 
This note is probably best read in conjunction 
with S. Giisken's plenary talk writeup sum- 
marising SESAM's attempt to identify sea-quark 
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Error analysis. One of the first of SESAM's 
results [D was to find that, given long enough 
trajectory lengths in the Hybrid Monte Carlo 
(HMC), we can obtain clear signals in the auto- 
correlation functions for simple observables (Wil- 
son loops, correlators). We find autocorrelation 
times which are much lower, i.e.O(25) trajectory 
lengths, than some of the conjectures to be found 
in the early HMC papers. The statistical errors 
presented in this paper are obtained from an anal- 
ysis of 200 configurations (per sea-quark) picked 
from a total trajectory length in the HMC of 
5000. A blocking analysis is performed and the 
quoted errors are always from a blocksize of 6, 
where we find errors to run into plateaus. I be- 
lieve this is the first time trajectory lengths in a 
HMC are long enough to present a reliable error 
estimate! Obviously, this is a vital step in any 
attempt to find signals of unquenching, particu- 
larly in the area of spectrum calculations since 
quenched spectra are (i) in agreement with ex- 
periment within 0(10%) (see e. g. ^) and (ii) 
can be measured to ever higher precision. 
Analysis with two degenerate light quarks. 
Next, I would like to discuss some issues raised 
recently by SESAM concerning the analysis of dy- 
namical configurations with 2 degenerate light 
Wilson fermions Q . The new feature of our anal- 
ysis is the use of spectrum data with unequal 
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valence and sea quark masses, in particular to 
determine quantities with strange quark content. 
We find our data to be well parameterized by a 
linear ansatz both in the sea and valence quark 
masses (see figure 5 in We have applied this 
method to the determination of the strange quark 
mass Q , which, so far ||] , has only been obtained 
in a sea of strange quarks. The effect is a sig- 
nificant increase in the two-flavour strange quark 
mass, whereas our dynamical light quark masses 
are much lower than the quenched ones. A sur- 
prising feature in the analysis is obtained when we 
attempt to study the light quark masses at fixed 
sea quark mass (this has been termed "partial 
quenching"). We find that valence quark masses 
need to be tuned to negative values to make the 
pseudoscalar masses vanish. In other words: the 
critical kappa-values at fixed sea quark lie below 
the true critical kappa of pions with equal strange 
and valence quark content. We have noted Q 
that light quark masses measured with respect to 
the partially quenched critical kappa values are 
much higher, practically in agreement with the 
values found in the quenched theory. 
Spectrum Results. The spectrum and decay 
constant results are displayed in figures I and 
H. We compare the dynamical results to (i) a 
quenched simulation we have performed at /3 = 
6.0; the two simulations are at comparable lattice 
spacings. To mimic the situation of the dynami- 
cal simulation, the quenched chiral extrapolations 
were performed linearly; (ii) the result one ob- 
tains (partially quenched) by working at a fixed 
sea-quark mass of k = 0.1575. Studying observ- 
ables with only light quark content, we see some 
improvement for the A, the decay constants and 
the J parameter, albeit within large errors. The 
nucleon, however, seems to be unimpressed by the 
presence of sea quarks. In the strange quark sec- 
tor the familiar problem from quenched QCD is 
still present; the K and the K* cannot both be 
matched simultaneously to experiment. 

In figure ^ I show the lattice spacings obtained 
from four different observables in light and heavy 
quark physics. Much better agreement seems to 
be obtained in the full theory. 
Finite volume errors. It was pointed out last 
year by S. Gottlieb M that little is known so far 
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Figure 1. Spectrum results. 
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Figure 2. Decay constants and J-parameter. 
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Figure 3. Lattice spacings from different observ- 
ables {F denotes the value obtained from the 
static potential at ro; the energy level splittings 
are obtained from lattice NRQCD |^.). 
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Figure 4. 
data. 



Finite volume effects in the spectrum 



about the size of finite volume errors in dynamical 
Wilson fermion simulations. Figure ^ shows the 
spectrum data from the SESAM and T^L collabo- 
rations at TO7r/wp = 0.69 obtained on lattices of 
16^ and 24^. While the errors are yet too large for 
definitive statements there is a general downward 
trend of the order of 2, 3 and 5% respectively for 
pseudoscalar, vector and nucleon masses. 
Quadratic contributions. Given that we have 
been using only three values of the sea-quark mass 
in our present analysis we have been restricted to 
linear chiral extrapolations. Particularly for the 
vector but also for the nucleon, the data seem to 
exhibit some downward curvature towards lighter 
sea-quark masses. Even with only three sea-quark 
masses at our disposition there is a nice con- 
sistency check for the curvature of the nucleon 
due to the Feynman-Hellman theorem which com- 
bines the spectrum analysis with our calculation 
of the TT- nucleon sigma term [FIS] . Given that 



iq {P\{uu + dd)\P) = m, 
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we can use our results for a to constrain the linear 
and quadratic terms of the proton mass. The 
result of this test is shown in figure ^ where the 
linear fit to the data is also plotted. The nucleon 
data are very well described by eq. |l|. The effect 
is to lower the mass of the nucleon in the chiral 
limit by 0(10%). It is a very pleasing result to 
find the spectrum data match up so well with the 
flavour-singlet operator calculation. 



Figure 5. Linear fit of the nucleon data and 
quadratic fit constrained by the a term. 



3. Conclusions 

We are currently in the middle of the analysis 
of the spectrum data of SESAM 's and TxL 's large 
statistics simulation of QCD with 2 dynamical 
flavours. We now have good control over statis- 
tical errors and have developed a new method to 
analyse particles with strange quark content in 
a sea of light quarks. We are pushing to lighter 
quark masses, mj^/mp « 0.5 and to larger lat- 
tices, enabling a badly needed finite volume study. 
Whereas the quark masses are found to be rather 
sensitive to the inclusion of sea-quarks a clear sign 
of unquenching is still missing in the spectrum 
data. 
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